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The problem of surface effects at a fluid boundary created by the force field of finite value is investigated. A 
classical simple fluid with a locally introduced field imitating a permeable solid is considered. The cases of micro- 
and macroscopically smooth boundary are examined and the analysis of static membranes is performed within the 
framework of general consideration. 

Henry constant of adsorption and its connection with Henry constant of absorption, specific surface Q- potential 
(grand potential) 7 and the surface number density are determined. High-temperature expansions for basic values are 
obtained. In the low-temperature limit it is shown that the results coincide with the previously considered problem 
of impermeable wall. 

"The surface tension coefficient" decomposes into a value proportional to Henry constant of adsorption depending 
on the position of the separating surface, and a universal nonlinear surface coefficient. 

Two approaches to this problem are analyzed: through the surface cluster expansion and through the pressure 
tensor. 

Within the first approach, the series in powers of activity is obtained for the surface part of the omega- potential. 
This expression is similar to the cluster expansion for pressure but in contrast to this case the integrals of Ursell 
factors contain the multipliers depending on the potential of particles interaction with the external field. 

Within the second approach, Kirkwood-Buff formula for 7 is extended for the case of the systems under considera- 
tion (the field of finite value). As a function of activity, depending on the situation, the surface terms may begin with 
either a linear term or a quadratic one, which corresponds to the presence or absence of adsorption in the classical 
understanding. 

It is demonstrated that the derivative of the tangential component of pressure tensor with respect to the chemical 
potential coincides with the near-surface number density (both averaged over the transition region), which, first of 
all, proves a complete identity of the approaches of "cluster expansion" and "pressure tensor" within the limits of 
their domain of existence, and second, gives the near-surface virial expansion which determines the exact equation 
of state of the "two-dimensional" gas of the near- surface region. 

The procedure of the "mechanical definition" of 7 is upgraded and substantiated. Coincidence of pressure acting 
on a transverse wall and the tangential component of pressure tensor, both averages over the transition layer as well 
as the symmetry of the solution with respect to the permutation of sorbent - fluid are demonstrated. 



1. Introduction 

Surface effects at the boundary between a fluid and a 
solid body considered using the methods of statistical me- 
chanics have been a subject of permanent attention of 
researchers for more than 60 years [T l7 
works are cited in [sllll. 



Some earlier 



Two approaches to this problem have been formed his- 
torically. The first one is based on constructing the sur- 
face cluster expansion - a series in powers of the activity 
for surface values [lj|5||8j[9] . The second one appeals to 
the notion of pressure tensor 



12-14,17 18 



In those works, a solid body is modeled by a constant 
field and the structure of the field is such that is has an 
impenetrable core. 

The reference points of this sphere are "surface tension 
coefficient" (the surface part of the omega potential or the 
grand potential) 7 and "gas adsorption on a solid surface" 



(including Henry constant of adsorption). 

Unlike in the case of hard core potential, investigations 
of a permeable wall are less thorough. As far as we know, 
works dealing with the statistical mechanical investigation 
of this problem had been completely absent until the 80- 
ies of the past century. The works that appeared in this 
area later use mainly the DFT approach (for example, 
see [19]). 

In our previous work 



20 we developed a consistent 
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approach to the problem of the interface of impermeable 
wall/fluid that allowed us to obtain closed expressions for 
key values. In addition, the equivalence of "tension"- and 
"adsorption" -based approaches was demonstrated. 

The present work spreads the approach of [20] for the 
case of a permeable wall. The work gives the general cal- 
culation of the omega potential (section |3| , specific sur- 
face omega potential 7 (sections |4j [7| and the surface 
number density (sections [S] and [?]). The upgraded proce- 
dure of the "mechanical definition" of 7 is given (section 



10 1 and Henry constant of adsorption is obtained (section 
6|T The equivalence of the cluster expansion and the ap- 
proach based on a pressure tensor is proved (section 13 1. 



Finally, the surface virial expansion is given in section [18| 
In addition, a connection between Henry constants of 
adsorption and absorption is established (section |6]) and 
the high-temperature expansion of surface values is pre- 
sented (section 15 1. Static membranes are considered in 
section 



17 and the symmetry of the problem in section 



2. Primary definitions and relations 

2.1. Canonical ensemble 

The probability density to find a given spatial configura- 
tion of a specific set of particles [ij p. 181] is given by the 
expression 
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Here N is the number of particles in the system, 
/3 = l/ksT, ks is Boltzmann constant, T is tempera- 
ture, t/f jY is the energy of the interaction of particles 
with each other, V is the volume of the system. Integra- 
tion is performed over the coordinates of the particles of 
the ensemble r}^^i...rjsi ■ Zpf is the configuration integral 



(2) 



Passing to distribution functions for an arbitrary set of 
particles, we obtain: 



(fc) 



(3) 



(k) 

where Qfj \ sets the probability density to find a given 
configuration of k arbitrary particles for the canonical en- 
semble. 

2.2. Grand canonical ensemble (GCE) 

Let us average equation ([s]) over the fluctuations of 
the number of particles, that is, let us apply operation 
Sw=o to both of its sides, where 
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is the probability for the GCE to have a definite number 
of particles N inside volume V. Here z is activity: 



(5) 



where fi is the chemical potential, A = /i/v^STmifesT, h is 
Planck constant, m is the mass of a particle, and Sy is 
the grand partition function of the system having volume 
V: 
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We obtain: 
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X / exp(-/3C/j^+^^Jdrfc+i...drfc+Ar|. 



1 fc specifies the probability density to find certain 
configuration of k arbitrary particles for the GCE. For 
ideal gas ^ = g'^, where g = N /V is the number 

density. 

2.3. Ursell factors and factors of partial localization 

(k) 

Ursell factors \. are also called cluster functions, it is 
these functions that appear in the known expansion of 
pressure in powers of the activity (41 p. 129], f2Tl p. 232] 

P(z, T) = zksT + kBTY,^ U^lkdr2-dru. (9) 

k=2 ■ 

They are also included into the corresponding expan- 
sion of the number density 

g{z) ^z + zf^^ f Uti%dr2...dr^+u (10) 

n=l 



which is evident from the relation Q 
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Ursell factors possess the locality property which is es- 
sential for this consideration: they decay rapidly with the 
removal of any group of particles including the unitary 
one. 

Taking into account that (|9]) is in fact the expansion in 
Taylor series, and keeping in mind the local character of 
Ursell factors we may write down the equation 



gkp 



(12) 



n=2 



which is true on the macroscopic scale and where the 
derivatives are taken at the point z = 0. Here S{r) is 
the Dirac delta function. We naturally assume that index 
k does not reach the macroscopic values. 

Some other properties of Ursell factors are described in 
Appendix |A.1| 

The factors of partial localization B^^'^^^j^ are the hy- 
brids of Botlzmann and Ursell factors. They are included 



^To avoid bulky equations, we will not indicate constant temper- 
ature; this will be implied for all the derivatives. In addition, 
differentiation parameters will not be indicated in the cases when 
their values are evident, for example are defined by the opposite 
side of an equation. 
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in the expansion of the distribution functions in powers 
of the activity 
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where 



1 (nev^) 
(n^l/) 



(14) 



are characteristic functions. 

Some properties of B^'l^^f, are described in Appendix 

EH 



2.4. Ursell functions 

UrscU functions are analogous to UrseU factors but they 
are composed not on the basis of Bohzmann factors but 



on the basis of distribution functions of the type ([8|) 22 
They are also called localized distribution functions. 
A definition of Ursell functions is the equality similar 



to ( A.l I and true for arbitrary k > 1: 



2.5. Presence of an external field 

The configuration integral of a nonuniform closed system 
is given by expression 



N 

zYj = j exp(-/3^t;, - PU^ j,)dr^...drN, 
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where Vi is the energy (= potential) of the interaction of 
the i-th particle with the field. 

For the GCE, wc introduce the quantity 



1 + E 
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which is evidently the grand partition function of the sys- 
tem in the presence of an external field. 

Using the properties of the fractional generating func- 



(101 



tion (|8| (see Appendix A. 2 B.ll we obtain an analog of 



oo „ „ n+1 

(20) 

where g{ri , z) is the number density in the presence of an 
external field, and 



/cq/ - — - k ^ 
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where {n} designates some partitioning of the given set of 
k particles with coordinates ri,...r/c into disjoint groups 
{tIq,}, I is the number of groups of a specific partition- 
ing, fcfj is the size of the group numbered a, summing is 
performed over all possible partitionings. 

These functions also possess the property of locality: 
they decay rapidly with the removal of any group of par- 
ticles including a single particle. Several first functions 
for the uniform medium look like 



(1) 



Q 



(16) 



(2) _ A2) 
1,2 ^ <Jl,2 



7r(3) 
1,2,3 



^1,2,3 



similarly to (A. 2 ) 



We will need the expansion for the uniform medium 



^ n 
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connecting them with Ursell factors, as obtained in [23] . 
It should be stressed that integration is carried out in 



( 17) over the whole space 24 . Note that another kind of 



Ursell factors, an equivalent one, was also used in 23 



= exp(-/3t)i) = exp[-/3w(ri)]. 



(21) 



and the analog of ( 13 1 



Q^rLio) = [ n + E It (22) 

/7n-\-k 
[ Jl di Sj";jj'_^j,(ir„j+i...dr„,+fe|. 

i—m-\-l 

3. Basic equations 

We will consider an open statistical system in the field, 
with the scope of the field smaller than the size of the 
system itself, and the position of the field at rather large 
distance from the boundaries of the system, that is, we 
talk about a situation more likely denoted as "the field in 
the system" . 

We will assume that the potential of the field is con- 
stant, except for the boundary regions, unless otherwise 
stated. Its changes on the macroscopic scale can be con- 
sidered admissible. 

The potential of external fields will be considered to 
increase/decay sufficiently rapidly to provide the conver- 
gence of the corresponding integrals when moving into 
the depth of the body or outside, while the interparticle 
potential decays sufficiently rapidly to ensure the conver- 
gence of the zero and first moments of Ursell factors (so, 
we are not considering Coulomb potentials in this work). 

In some cases we will accept that the region affected by 
the field of the solid body has a finite radius, and speak 
of the transition region between the solid body and the 
fiuid. This approximation in most cases does not have a 



3 



principal character and will be used to simplify consider- 
ation. 

Finally, unless otherwise stated, we will accept that the 
region of the transition from the solid body to the fluid is 
much smaller than the linear size of the field, and the pro- 
file of the field along its gradient is the same everywhere. 

Some sections of this paper are similar to the corre- 
sponding sections of 
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Taking the logarithm of the ratio of the partition func- 
tions (19 1 and Q, we obtain 
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where we have used a logarithmic type of the generating 
function of the Ursell factors. 

In the above-indicated assumptions, we may remove 
symbol V from the integrals over the coordinates of parti- 



cles in ( 23 1 and accept that integration is carried out over 



the infinite space. 



Equation (23 1 is transformed as follows: 
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where 

z' = zejq>{-l3vo) (25) 

is activity, and vq is the potential of the external field 
inside the region affected by the field at the sufficient re- 
moteness from the boundaries, and where designations 



and 



exp {-I3vi) - exp {-I3vq) 
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1 - exp (-/3wo) 
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(27) 



are introduced, generalizing Boltzmann factors (21) and 
Mayer functions 



LPi exp (-/3wi) 

respectively. One can see that 



(28) 



(29) 



and functions 9i and ipi decay inside and outside the field, 
respectively, and tend to 1 in the contrary case. 
In addition, they have the property 



+ z ipi = z6 



(30) 



(Symmetry property of these functions (120), see sec- 
tion 
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The first sum in right-hand side of ( 24 1 , as one can see 
from ([9|, is proportional to the difference of pressures of 
uniform systems. 



The second term in the right-hand side of (24) is pro 



portional to the area of the surface separating the field 
region from the rest part of the system. Indeed, its struc- 
ture is such that at least two particles always stay on 
different sides of the "field boundary" . If we expand the 



product in the arbitrary term in (24 1, we will obtain the 
expressions of the type of 



(31) 



where l<A:<i — 1; 1 < I < t — 1; k + I — t. Assume 
that the first particle belongs to group 9i. (This may be 
made without loss of generality due to the symmetry of 
Ursell factors under permutations). With the fixed first 
particle, due to the locality of Ursell factors, integrals of 
the type of ( [sT] ) are defined by the local area near it. 
As a consequence, / . • . dr2 ■ ■ ■ drt are independent of the 
displacement of the first particle along the boundary of 
the system. 

With the displacement of the first particle from the 
boundary to the depth of the field region, they decay 
rapidly due to factor 9i. With the displacement of the 
first particle from the boundary to the depth of the free 
fluid, they decay rapidly due to fixing factors (pj and the 
local nature of Ursell factors. Integrating along the sur- 
face, we can factor out of the sum the area and make other 
simple transformations to obtain 



-P{z)V +[P{z) ~ P{z')] / ip{r)dr + i^{z,z')A, 
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(32) 

where ^l^ is the omega- potential of the system with the 
field immersed into it, A is the area limiting the field 
region and where we introduced designation 



OO 

i^{z, z) = ksT^- 
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Oiz'' + - Y[{ze, + z'^i) U[';>^dxidr2...dr, 



Here is the coordinate directed along the field gradi- 
ent, and the direction of xi axis is chosen so that dxi > 0. 

To remind, integration is carried out over the infinite 
space. 



In (33), integration over xi may be performed. Chang- 
ing the variables x'l = xi, r[ = Ti — ri, i = 2, . . . k and 
using the fact that the Jakobian of this transformation is 
equal to 1 and the invariance of under translations. 
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we obtain 



oo 

v{z, z') ^ksTY,-^ 9i'::kM':Lkdr2-drk, (34) 
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where 
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Let us denote the volume of the constant field as Vi, 
and the volume of the transition region as Vt ■ As we have 
already mentioned, Vt <^Vi. 



At the first glance, from ( 36 1 we obtain three versions 



of accounting for surface effects for different values x': 

n = -Piz){v ~Vi- Vt) - Piz')iVi + Vt) 

~ [Piz) - P{z')]A [ 9{x)dx + i^{z,z')A, (37) 



v) + z'(p{x)] Y\[zO{x + Xi) + z'lp{x + Xi)]^dx 



is the function of variables X2, ---Xk, symmetric under per- 

(k) 

mutations. Evidently, unlike U-^ \ it is not invariant un- 
der translations because it depends on the localization of 
the external field. 



It should be noted that the variables in (34) are sep- 
arated: the external potential is present only in g^^~^\ 
while the interparticle potential is present only in U''^\ 

We will see below (section [l4| that coefficient v is sym- 
metrical under permutation of z and z', with the corre- 
sponding changes of the level of energy readout and its 
sign. 



Equations ( 32 1 and ( 34 ) generalize expressions for the 

to the 
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f2- potential and nonlinear surface coefficient 
case of the finite field. They are basic for subsequent con 
sideration. The first term in the right side of ( 32 1 has 
purely volume properties, and the third one has purely 
surface properties. The second term as is easily seen from 
its structure possesses both volume and surface proper- 
ties. 

As can be clearly seen, the expansion of a function 
v{z, z') in series of z, z' begins with the quadratic form. 
Other representations of this quantity will be established 
below. 

4. Variants of accounting for surface effects 



Equation (32) may be transformed to the form: 



VI = -P{z){V -V) - P{z')V' 

x' 

- [P{z) - P{z')] \ j 9{x)dx - I ^{x)dx 



A 



+ v{z,z')A, (36) 

where x' is some arbitrary point inside the transition layer 
or near it, determining the volume V , and x is the coordi- 
nate directed along the field gradient. (Equally well, the 
regions of integration in equation ( |36[ ) can be bounded by 
the limits of the transition layer because the integrands 
become zero outside it.) 

We see that there is a definite degree of arbitrariness 
in partitioning the volume and surface terms. In reality, 
this freedom exists only within the transition layer or near 
it. Indeed, one can easily understand that otherwise we 
speak of the trivial compensation of two identical terms 
having opposite signs. 



-P{z){V -Vi)-P(z')Vi 
+ [P{z) - P{z)]A I ip{x)dx + v{z, z)A, (38) 



n 



n = -P{z){V - Vo) - P{z')Vo + v{z, z')A, (39) 

where It is the length of the transition region, and Vq is 
the volume of the body limited by the surface that is set's 
by condition 



9{x)dx — / ip{x)dx 



(40) 



or 



[exp {—(iv) — exp {—(3vQ)\dx — [1 — exp {—(3v)]dx 



(41) 

In the left side of this equality, integration is carried 
out from the field region outwards, while in the right part 
it is made to the depth of the fluid, and it determines the 
position of the surface of interest for us xq. 



The variant ( 37 ) corresponds to x' being at the bound- 



ary of the uniform fluid, that is, relates the surface effects 



to the solid body (fleld), while ( 38 ) - to x' at the boundary 
of the uniform solid body, and relates the surface effects 



to the liquid. Variant (39) makes the linear surface terms 



equal to zero using the property ( 40 ) . Depending on the 



shape of potential, x' may be either inside the transition 
region or outside it. 

In reality, due to the symmetry of the problem of per- 
meable wall with respect to the permutation of the fluid 
and field (see section 14 ) , we have only two different ver- 



sions. The variants (37) and (38) are equivalent in this 



case. This can be easily proved by making substitutions 
in ^ a.sV{ ^V-Vi- Vt and (fill) - (fm I. 



Of course, other version based on equation ( 36 ) can 
be considered. The results obtained here we will discuss 
below. 



Expression ( 36 ) is equivalent to 



-P{z){V - V) - P{z')V' + jA, 



(42) 



where 



7 = -[P(z) -P(z')][ J 9{x)dx- J ifiixjdx'^ +iy{z,z') 

—00 x' 

(43) 

is the general form of the specific surface rj-potential. 
This expression is separated into the linear and non-linear 
(v) terms with respect to pressure. 
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5. Surface number density 



Differentiating ( 36 ) with respect to the chemical potential 
and taking into account the fact that |^ 



= - 



(44) 



V,A 



and keeping in mind ( 11 ), we obtain 



N = g{z){V - V) + g{z')V' 



+ [p(z) — g(z')] / 9{x)dx — I ip{x)dx 



A 



-^/3z^^.(45) 
oz 

Differentiating z') with respect to z it is evidently 
necessary to consider Z SiS Si function of z in accordance 



with (25) 



Expression ( 45 ) can be written in the form 

N^Nb,+Nb,+N„ 



(46) 



where A'^^i , and are the numbers of volume par- 
ticles in the different phases and the surface particles, 
respectively. For them, we have the expressions 

= g{z)iV - V), = giz')r, (4?) 



N, = [g{z) - g{z')] / e{x)dx - / ip{x)d 



A 



Apz^-^^.m 



Passing to the surface number density 



A ' 



we obtain 



Qs — \_q{z) — g{z')\ j 9{x)dx— / tp{x)dx ~l3z 



(49) 



di'{z, z') 



dz 



(50) 



Taking into account that from ( 20 ) follows that 

oo 

oz 



dx, 
(51) 



g{r) - 0{x)g{z) - ip(x)g{z') 
we obtain finally after simple transformations 

x' oo 

Qs = / [q{x) ^ Qi^')] dx + [g{x) — g{z)] dx, (52) 



where g{x) is the number density in the vicinity of the field 
boundaries depending on the coordinate directed along 
the field gradient. 



We will omit the averaging sign for A'^ of open systems. 



We will see further on (section [9| that ( 52 ) agrees with 
the usual expression 



dfi' 



(53) 



Formula (52 1 generalizes the expression for the surface 



number density 20 to the case of the field of finite value. 



Below we obtain the second form for gs ( 72 ) , in which 



the terms are grouped according to linearity criterion. 



Expression (52) has a simple sense. Relying on the 



evident relation 



iV = y g{r)dr 

V 



(54) 



and using (46) and (47), we obtain 



g{r)dr - g{z){V - V) - g{z')V' 



(55) 



[g{r) - g{z')] dr + / [g{r) - g{z)] dr, 



v-v 



which under the above assumptions is identical to ( 52 ) . 



So, this expression is true for any densities of the fluid. 



The correspondence between ( 52 ) and ( 55 ) provides evi- 



dence of the internal consistency of the approach applied. 

To conclude, it should be noted that, unlike in the case 
of the equilibrium system liquid/gas, similarly to the case 
of impermeable wall [20| , superimposition of the condition 
Ns — 0, as one can see from ( [52| , results in the depen- 
dence of the position of partitioning surface on number 
density. So, in this case, too, the physical sense of such 
a partition is absent. Instead, in some cases it is possible 



to use the condition ( 40 1 , coming to the variant of "zero" 
adsorption. 

6. Henry constants of adsorption and 
absorption 



For Henry constant of adsorption 



Kh = lim 



!?->-0 g[z) 



(56) 



from ( 50 1 or ( 52 ) we obtain 



X 

Kh — J [exp(— /^tj) — exp{—^V())]dx 



oo 

J [I^exp{~-I3v)]dx. (57) 



This formula generalizes the previously obtained ex- 
pression 



20 to the case of the field of finite value. 



The dependence of on x' is universal and can always 
be extracted in the explicit form. Indeed, subtracting one 
value of the constant from another we make sure that 

Kh{x') = (x' - x")[l - exp{~l3vo)] + Kh{x"). (58) 
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If we choose the value x" = Xq, determined by (40 1 or The relationship (65) corresponds to high temperatures 



(41 1, then, because 



ifff(xo) =0 



we obtain 



Kh{x') = \x' - xo(T)][l - exp(-/3«o)] 



(59) 



(60) 



and the problem of K^i determination is reduced to the 
calculation of xq. 

The sense of expression (60) is quite clear: this is a 



coefficient proportional to the difference of densities in the 
limit of low activities: [1 — exp(— /3wo)], multiplied by the 



volume of the parallelepiped with thickness 
and unit base area. 



It follows from (40) that 



6{x)dx + / ip{x)dx 



(61) 



The surface passing through point xq by virtue of ( 59 ) 
can be called the surface of zero adsorption. 



For Henry constant of absorption, using expansion ( 20 ) 
we obtain 



fcff = lim 



g{z') 



e-i-O g{z) 



exp(-/3uo)- 



(62) 



Equations (57), (62 1 determine, within the framework 



of the approach applied, the connection between the two 
Henry constants: adsorption and absorption. Here we 
speak of the parametric dependence on temperature under 
the fixed potential. 



Indeed, excluding temperature from equations (57), 



( 62 ) we obtain 



K 



H 



dx. 



where 



v{x) 



v{x) 



(63) 



(64) 



is the potential of particle interaction with the external 
field, normalized so that it is equal to 1 inside the body. 

The procedure used here may seem somewhat artificial, 
but the key factor is that the reduced potential v in fact 
does not contain parameter vg as a consequence of binding 
the actual potential at the boundary of the body to the vq 
value. In other words, parameters vq and /3 are included 
in expressions only in the combination f3vo and thus the 
procedure is legal. 

Now we can make some general conclusions concerning 
the mutual behavior of and kjj. The relations that 
we have in mind are of the form 



Kh 
Kh 





x' 



X2. 



(65) 
(66) 



As one can easily make sure, these relations follow di- 



rectly from ( 63 ) 



and agrees with the simple fact that the surface effects 
disappear in the uniform medium. 



The relation (66) is true, provided that v{x) > in the 
whole space. This means also that > 0. The value 
X2 corresponds to the point in which the external field 
becomes zero. So, this expression is true only under the 
assumption of the finite thickness of the transition layer. 
(Due to the symmetry with respect to the substitution 



field O fluid (section 14 ) the case under consideration in 
fact includes also the case v{x) < 0. That is, here we 
exclude only the variant of sign change for v{x).) 



The expression (66) has a simple sense: for x' > X2 



this is the volume of a parallelepiped that has the base of 
the unit area and the height x' — X2- Multiplying it by 
the volume number density we obtain the surface number 
density. (Of course, we consider the limit of low densities.) 
The sense of Kh for x' < X2 is clear too: this is the 
negative amount compensating the fictive amount of fluid 
in the region of the solid body. 

So, provided temperature tends to zero in the point 
where the field vanishes (X2), in fact, an infinite wall ap- 
pears. 

It is clear from the above considerations that the equa- 



tions ( 66 ) give the maximal value of constant Kh for non- 



negative potentials under x' > X2. 



Thus, expressions (65) and (66) show that adsorption 



and absorption anticorrelate in the case of the absence of 
the regions of negative potential. 



In (63), terms with x' may be distinguished; then 



00 

Kh = x'{l -kH)+ J [fcf - kH]dx - j kl^"^]dx. 
-00 

(67) 

As an example, consider the case of the linear change 
of potential in the transition region 

{1 {x<xi) 
{X2 ~ x) / {X2 — Xi) (xi < X < X2) , (68) 
{X2 < X) 

where Xi and X2 are the coordinates of the beginning and 
end of the transition region from the body to the fluid. 
Carrying out an elementary integration we obtain 

1 ^ kn 

Kh = x'{l - kn) - — (2^2 - a^i) + xikn - X2. (69) 
m kn 

It should be stressed once more that the above consid- 
eration relates to the parametric dependence of Kh on 
kn- For the direct dependence through vq, the situation 
is more complicated because it is not clear how to model 
the behavior of the potential v{x) in the transition region 
under variations of vq. From the viewpoint of expression 



(60), the difficulty is that xo depends on vq in a non- 



trivial manner. In all likelihood, this dependence can be 
considered only within the frameworks of specific models 
the significance of which is unobvious. 
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7. The basic form 



8. Macroscopically smooth field 



It follows from (36) that 

n= -p{z){v-v')~-p{z')v' 

- [P(z) - P(z')] {x' - xo)A + iy{z, z')A, (70) 
or, by virtue of ( [42| ) 

7 = - [P{z) - P(z')] {x' - xo) + y{z, z'). (71) 

This equality is one of the basic expressions for the con- 
sideration of surface phenomena in this case. It binds four 
main surface parameters: specific surface f2-potential 7, 
the position of the planes of zero adsorption xo and sepa- 
rating surface x' , nonlinear surface coefficient v and vol- 
ume parameters - pressure at both sides of the interface. 
This equality also generalizes the analogous equality 20 
to the case of the field of finite value. 



Let the characteristic distance at which the field under- 
goes essential changes be macroscopic. Substituting ( 12 ) 



into ( 34 ) and integrating over r2, ■ ■ ■ ,rk, we obtain 



Kz,T) 



^ z'' d'^P 

k=2 



|6l(a;) + <^(a;) exp(-/c/3wo) - [0ix)]''^dx (78) 



i^iz,z')= / 9{x)Piz) + ;p{x)Piz') - P{9z) dx. (79) 



Differentiating ([Tlj) with respect to the chemical poten- Substituting this expression into similarly to the 



tial and taking into account (51 1 we obtain 



case of infinite field 20 , we arrive at a logical result 



Qs = [q{z) - q{z')\{x' - Xo) 



+ 00 



[q{x) — 9{x)q{z) — ^{x)Q{z')\dx (72) 



- one more expression for the surface number density in 
the case of the arbitrary position of the dividing surface. 
Here two parts are clearly separated: one of them is linear 
with respect to density and the other is nonlinear. (This 
statement is true if we consider z' — z' {z), and the exter- 
nal potential is fixed. If we accept z and z' as independent 
variables, then the first term of the right part of ( [72| ) is 
evidently a linear form with respect to density, while the 
second term starts with the quadratic form of q{z), q{z'), 



= - J P{9z)dr ^- J P{ze-^'')dr 

V V 

P[li - v{r)]dr = - j P{r)dr, (80) 

V V 

where we used a known expression for the chemical po- 
tential in the force field. 



We obtain from (61 ), ( 71 ) and (79 1 



+ 00 



7(a;')= / [P{z')-P{x)]dx+ / [P{z)- P{x)]dx. (81) 



which follows from (51 1 and (34), (35).) One can easily 



make sure that ( 72 ) coincides with ( 52 ) 



So, we may consider the separation of surface values 
according to the linearity principle on the same ground 

as that for the localization principle. Indeed, if we define 9. "Local pressure" 

the linear surface number density as 

9sd = [g{z) - g{z')]{x' - xq), 

and the nonlinear one as 



A characteristic structure that appeared in (81 1, as will 
be seen below, is universal for 7. 



(73) 



+00 



Integrating ( 44 ) over the chemical potential we obtain 

n^- I Ndn', (82) 



[g{x) - 9{x)g{z) - ip{x)g{z')] dx, (74) 



analogs of ( 53 ) 



then, taking into account (11) and (51), we obtain the 



where we omitted an arbitrary function of volume, tem- 
perature and field appearing as a result of integration be- 
cause the f2-potential should tend to zero in the low den- 



d\[Piz)-P{z')]ix'-x„)} 

gs,i = ^ ^ (75) 



sity limit. Using expression (54 1, we come to the relation 



and 



gs,n 



9/i 

di'{z, z') 
(9/i 



(76) where 



according to (531, (71). To remind, in (75 1, (76) we dif- 



ferentiate both 2, and z' . Naturally, 

gs = gs,i + gs,n 



^ = - j P*{r)dr, 

V 

P*(r)= J gir)dfi', 



(77) 



and expressions ( 73 ) - ( 77 ) determine the surface number 



density in the most general form. 



which is equivalent in this case to 

dP*{r) 



dfj, 



(83) 



(84) 



(85) 
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This means that the number density for the system in 
the field and P* are connected through exactly the same 
relation as usual density and pressure for the uniform sys- 



tem (111. 



The value equivalent to P* was considered for the im- 



20 and we will also call it "local leads to this definition 25 



permeable wall 16 
pressure" . 

Under the condition of Mayer- type cluster expansion 



10. "Mechanical definition" of 7 

The procedure of compression and extension of the sys- 
tem in two different directions with overall conservation of 
the volume and with changes of definite internal surfaces 

p. 43]. This procedure will be 



validity, we obtain from ( 20 ) 
P*(ri,z) = ksTzOi 



(86) 



n=2 ■ •' i=2 

which also resembles the standard expansion of pressure 
in powers of the activity Q. 

One can see from (86) that P* rapidly tends to P as 



applied in the upgraded form to the field/fluid interface. 
The consideration given below is a generalization of the 
procedure used for the case of infinite field [20j . 

In virtue of fundamental identity of the boundaries, it 
is necessary to take into account surface effects not only 
on the interface that we are interested in, but also at the 
boundaries of the system. 



I 



the distance from the field sources increases. 



It should be stressed that equation (83 1 is absolutely 
general and true for the system in the force field of the 
arbitrary configuration. In the case of macroscopically 



smooth fields it turns into expression ( 80 ) 



Similarly to the above consideration for the density of 



particles, equation ( 83 1 can be brought by means of simple 



transformations to the form 

n= -p{z){v -v')-p{z')v' 



(87) 



J [P{z')-P*{r)]dr+ J [P{z) - P*{r)]dr. 



v v-v 
Integrating over the surface we obtain 

n= - Piz){v -V')- Piz')v' 

x' 00 

+ A [ [P{z')-P*{x)]dx + A [ [P{z) 



(88) 
P*{x)]dx. 



As it should be, the derivative of (88 1 with respect to 
the chemical potential gives (up to a multipliers) the ex- 
pression for the volume number density ( 47 1 and the sur- 



face one ( 52 ) 



For specific surface il-potential 7 we obtain 

x' +00 
— 00 x' 

(89) 

Equation (|89|) generalizes the expression obtained pre- 



viously 16 , 20] to the case of the finite value potentials. 



Certainly, expression (42) remains true in this case. 



Differentiating ( 89 ) with respect to the chemical poten- 



tial and taking into account ( 52 ) , we make sure that ( 53 ) 
is also true. 



From (33 1 taking into account ^ and (86) one can 
easily obtain an expression for the nonlinear surface coef- 
ficient 



iy{z,z')= J \^e{x)P{z) + if{x)P{z')- P*{x^dx, (90) 

— CO 

determining the form of v in terms of "local pressure". 



The similarity of expressions (90) and (79) should be 
stressed. 



hi 





1 























t 



Figure 1: The system for the mechanical determination of 
the surface il-potential of a fluid in a force field. 
The interface plane is singled out. The potential 
energies of the particles of the upper (index 1) 
and lower (index 2) parts are different. 

We will consider the limiting physical surface to be in- 
side the system at a microscopically large distance from 
its boundaries and at the same time this distance is much 
smaller than the size of the system. 

The procedure involves the following (see Fig[l]). At 
first we move the upper plane of the system downward by 
dhi and the lower one by 6h2 upward. (These quantities 
are, obviously, the absolute values of displacements.) The 
work done by an external force on the system is 



SWi 



h 

P^{l-2kf+Al J Pstjx 


It 



6hi 



P2{l~2hr+Al J Pst.dx 




5h2, (91) 



where Pgt. is the pressure in the region of the narrow 
strips near the edges of the system |^ differing from pres- 
sure at the center (in the bulk). Index s reminds that, 
generally speaking, this pressure does not coincide with 



^For simplicity, we accept that all the lengths of the transition 
regions; lt,lst,l* ,lv, connected respectively with Pt, Pst, P* ,0 
are the same. One may consider on default for averaging inte- 
grals that the integration is performed over the maximal of these 
lengths. 
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the tangential component of pressure tensor 20 

Pst + Pi- 



rn 



At the second stage, we ahow the system to expand but 
in a lateral direction 

5W2 = -[Pi{l-'2lt){hi~lt) 
h 

+{1 + 2/ii) / Pst^dx + P^il - 2k){h2 - 2lt) (93) 



h 



h 



{l + 2h2) J Pst,dx + l J Pstdx 





61. 



Here Pst is the tangential pressure along the boundary 
between the phases. The transition region is completely 
included into volume No. 2 (we speak here of the point 
separating heights hi and /12). 

Similarly to the case of infinite field, the condition for 
the conservation of volume is not quite trivial: we must 
consider the constancy of the volumes limited by the 
planes of arbitrary localization, set by coordinates x for 
system boundaries and x' for the phase interface {x' will 
be directed upward from the origin of the transition region 
in the second phase). It is clear that each volume should 
be conserved constant separately in order to compensate 
the volume terms in the final expressions. 

{I - 2xfShi = {l- 2x){hi -x + k- x')5l 

{I ~ 2xf5h2 = {l- 2x){h2 -x-lt+ x')5l (94) 



In view of the assumptions made, the lower and upper 
integration limits may be considered to be equal to —00 
and +00. So, it may be concluded that the expression for 
elementary work of change of the boundary surface area 
between the phases is 

x' +00 

5W = 5a\^ J [P{z')- Pst]dx+ J [Piz)-Pst]dxy 

— 00 x' 

(98) 

Because the process takes place at constant A'' and T, 
this work should be equal to the change of the free energy 
F — F{V,T, N). An apparent complication during the 
process is pressure variation; generally speaking, pressure 
changes as a consequence of adsorption. Varying 

F = ~P{z) (V - V) - P{z')V' + 7A + ^i{Nb, + ATfc, + Ns) 

(99) 

and taking into account the fact that 5Ni,^ +SNi,,, = —SNg 
and SV = SV = , we obtain 

SF = -{V -V')6P{z) -V'SPiz') (100) 
+ 7M + AS-f + {Nb^ + Nb,_ + Ns)Sn. 

In the right side, all the terms are cancelled pairwise 
except the third term, so finally 



6F = 7M, 



(101) 



as it must be if 



+00 



7 = y [P{z')-Pst{x)]dx+ J [P{z)~Pstix)]dx. (102) 

— 00 x' 

So, the "mechanical definition" of 7 also comes to the 



or conserving the terms of the first order of smallness, structure that we already know (81 1, (89) 



TT /7 ,61 

dhi = [hi H — —X — X + It — X 

2/i2 , ,.6l 

dh2 = (/i2 + —j-x -x-lt+x )y. 

Then, in the same approximation, we obtain 
SW = SWi + SW2 



(95) 



11. Tangential force and "local pressure" 

Now we will demonstrate the consistency of approaches 



(|89|) and (|102|). We differentiate (|99|) with respect to the 

,dP{z') 



area, where we use for 7(89) 

'dF 
dA 



Ml 



+6A2 



N,V 



Pst^dx+ {Pi-Pst,)dx 



+7 + A|l + (iV. 



dA 

+ Nb, 



V' 



dA 



Nb.)j^. (103) 



Pst,dx+ {P2-Pst,)dx (96) 



One can see that the first two terms in the right side 



SA 



{P2 - Pst)dx + {Pi- Pst)dx 



of (1031 are cancelled with (A"b^ + Nb2)dii/dA, and the 
fourth one with Ngd^/dA. We can easily make sure of 
this by passing under the integral to the derivative with 
respect to the chemical potential and taking into account 



where 



(|52|, so that finally we obtain 

[P{z')~P*{x)]dx+ / [P{z)-P*{x)]dx. 



5Ai 2hi)Sl; SA2 2h2)Sl 6A = ISl (97) 

are changes of the corresponding areas during the process. 
Similarly to the case of infinite field, the presence of free 
parameters (/ii,/i2) allows us further on to consider that 
only one surface is varied. 



OF 
dA 



+00 



N,V 



(104) 

So, we may identify ( 89 ) and ( 102 ), and thus we obtain 



P.fdx = P*dx 



(105) 
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or 



where the bar means averaging over the transition layer. 

Under the real conditions, in the absence of clear 
boundaries of the transition layer, the accuracy of equal- 
ity ( 106 ) is determined by the closeness of Pst and P* 



to P{z) and P{z') at the boundaries of the averaging in- 
terval. The nontrivial part of the equality is determined 
exactly by the difference of these values from their lim- 
iting values. So, the value of averaging interval in ( 106 ) 



must determined by the compromise between these two 
factors. 

12. Pressure tensor and near-surface 
distortions 

Now we will prove a similar equality for the tangential 
component of the pressure tensor 



P^P* 



(107) 



We will rely on the equivalence of the problems of in- 
homogeneity in the field and inhomogeneity in the region 

The initial expression will be 



of liquid/gas interface 11 



Ptixi) = giksT 



1 



(y2 ~ yi)Qi 



(2) dui: 



dy2 



dr-2 



(108) 



which is identical to that used in numerous works [3 12 



ITpIS . Here u^- is the energy of the interaction of i-th and 



j-th particles, y is the coordinate parallel to the surface, 
and the gradient of the external field is directed along x, 
similarly to the previous cases. 

A brief outline of proof is as follows. 



We integrate (86) and (108) over the coordinates of the 
first particle, then expand the subintegral expressions in 
the series of the activity using (22). Finally, equating 
the coefficients at the same powers of z, we prove their 
identity. 

The equality of the integrals of Pt and P* over volume 
is proved in Appendix [C| 

Now we will prove that Pt — P* — P va the region of 
homogeneity. Because all the directions are equivalent in 



this case, we may write down equations similar to ( 108 ) 



for X and z. Summing these three equations we obtain 



3Pt = igkuT — 

+ (2/2 - yi 



1 

2 

dui2 
dy2 



Q12 {X2-Xi) 
+ {Z2 - Zi) 



dui2 
dx2 
dui2 



dz2 



dry. 



(109) 



Collapsing the bracketed expression we obtain a well 
known expression for pressure 



Pt = 3kBT-\j B^^rJ^^dr2 = P 



p. 190]. 

It is evident at once for P* that for Oi — 1, i = 1, 2, . . . , 
([Sef is reduced to (|9|, while for Oi = exp (-/3wo) to P{z'), 
that is, 

P* = P{z) or P* = P{z') (111) 



in the region of homogeneity. 



(106) So, in view of the validity of equalities (110), (111) in 



the region of homogeneity and the equality of the integrals 
of Pt and P* over the system volume, we may conclude 



that ( 107 1 is true. 



Similarly to the case of Pst, the region of the transi- 



tion layer in ( 107 1 is determined by the closeness of 9 to 



exp (— /3wo) and 1 at the edges of the interval. 

There are some reasons to suppose that the local equal- 
ity is wrong, that is, generally speaking. 



Pt^P* 



(112) 



but a rigorous proof of (112 1 for the systems under con 



sideration goes out of the scope of the present work. 

We will not consider also the questions connected with 
the ambiguity of pressure tensor because we are interested 



only in its zero moment 17 



So, in view of (851 and (107), we estabhshed a con- 



nection of the tangential force (per unit length) with the 
number density near the surface, and thus with the be- 
havior of near-surface distortions 



d 



g(x)dx = — / Pt(x)dx (113) 
J 

It h 



dPt^ 



(114) 



where the correspondence to the macro-relation (11 1 also 
takes place. 

At the same time, we obtained an equality for the tan- 
gential force acting on the transverse wall 



P.f = Pt 



(115) 



which cannot be considered self-evident in advance. 



Finally, expression ( 89 1 can be also written in the form 



+00 



7(a;')= / [P{z')- Pt{x)\dx+ / [P{z) - Pt{x)\dx 



(116) 

This expression generalizes the Kirkwood-Buff equation 
for the arbitrary localization of the dividing surface [20j 
to the case of the force field of finite value. 



13. Surface cluster expansion and pressure 
tensor approach 



Equality ( 107[ ) allows us to prove the identity of cluster 
expansion and pressure tensor approaches in this general 
case too, similarly to the case of the impermeable wall 



20 . For this purpose, it is sufficient to demonstrate that 



(110) (43) with V in the form of (33) is equivalent to (116) 



Equalizing these two expressions and canceling identi- 
cal terms we obtain the condition of their identity in the 
form 



v{z,z') 



6{x)P(z) + ^{x)P{z')~ Pt{x) dx. (117) 



11 



Comparing (1171 with (90) and taking into account 



( 107 1 we ascertain that this equahty is vahd. 



Taking into account ( 71 ), we obtained one more expres- 



sion for the specific surface 57-potential 
^ = -[P{z) - P{z')]ix' - xo) 



Similarly, considering the surface number density (72) 
we obtain 

Qs{0,^, g{x)\z, z' ,x') ^ gs{0* ,(p* , g*{x*)\z* , z'* ,x'*), 

(125) 

where g*{x*) = g(x). 

Finally, Henry constants are to be considered. It follows 



6{x)P{z) + ip{x)P{z') - Pt{x) dx, (118) from the first two equations (119) that 



of course, coinciding with expression (116) but again ex- 



plicitly separating the terms linear and nonlinear with 
respect to pressure. Replacing Pt in (118) by P* we re- 



turn to the approach from the side of "adsorption" , while 
replacing it by Pgt we come back to the "mechanical def- 
inition" of 7. 

14. Symmetry of the problem 

Unlike for the problem of adsorption when we have an 
impermeable wall, the problem of absorption is principally 
symmetrical. Varying the parameters we may consider 
one or another region of the fluid be free or situated in 
the field. 

Indeed, let us make the substitution 



z = z 
z'* =z 

X — —X 

V* (x*) ~ v{x) 



(119) 



V, 







-Vq. 



This set is equivalent to the original physical context, 
but mirrored with respect to a; = and with a shift of the 
energy reference level. 



Performing the substitution of (119) into (26), (27), we 
see that 

0[v{x),Vo] = d[v*(x*) - VQ,-Vg] = ifi[v*{x*),VQ] 

if[v{x),vo] = !p[v*{x*) -Vo,-Vq] = 0[v*{x*),Vq]. (120) 



So, as one can see, for example, from (34), (|35|) coeffi 



cient v is invariant with respect to transformation (119) 



1/(0, Lp\z, z') = iy{ip* ,9*\z'* , z*) ~ 1^(6* ,(p*\z* , z'*), (121) 
where we have introduced the notations 



r = e[v*{x*),v^] 

If = ip[v (x ),va\, 
and the dependence of on 0, is functional. 



(In the derivation of ( 121 ), one must use the invariance 



of Ursell factors with respect to the change of coordinates 
sign.) 



In turn, it follows from (61) that 

* 

Xq — ~Xq, 



(123) 



which, on substitution into ( 71 ) together with other pa- 
rameters, gives 



^{9,^\z,z',x')^j{e*,^*\z*,z'*,x'*)- 



(124) 



g ^ g 
I* 

g = g. 



(126) 



(62) 



Then, for Henry constant of absorption we have from 

1 



h* — 



(127) 



while for Henry constant of adsorption from (56) or (57) 

K*h{x'*)z* = Kh{x')z. (128) 

Another property determined by the symmetry has al- 
ready been considered in section [4] 



Summing up, we may state that equations (121), (124) 



and (125) are control ones from the viewpoint of checking 



the problem on symmetry. 

For small densities, the surface number density is de- 



termined by gs^i (73 1. Demanding gsj > 0, we see that 



for g{z) > g{z') we must take x' > xq- 

So, if we accept that the denser medium is the fluid, it 
is the assignment of the transition layer to the solid body 
that agrees with the notion of adsorption as it is usually 
understood. 

15. High temperatures 



Using expression (61) and considering (3v{x) <^ 1 within 
the entire range of x, we obtain 



Xo = 



v(x) 



dx 



vlx] 



dx 



(129) 



So, we see that xq is temperature-independent at high 
temperatures. It follows for Henry constant of adsorption 



from (60) that 



Knix') = l3voix' - xo) + • • • - m^') + 



(130) 



(122) which corresponds to the limit (65) and where 



x')— / [vo — v{x)]dx — / v{x)dx. (131) 



The linear part of 7 is reduced to 

^ [P{z) - P{z')]ix' ~ xo) (132) 
= -g{z)vo{x' - Xq) -\ = -g{z)(l){x') + 

where we used the expansion of P{z') in Taylor series near 
point z. 
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Then, for linear with respect to the density part of the 16. Low temperatures 



surface number density ( 75 ) we have 

0s,i = g^vo{x' - xo)>ct -\ = g^>CT(t){x') + . . . , (133) 

where 



ldQ_ 
~Q~dP 



(134) 



is isothermal compressibility. 

Now the nonlinear part of 7 will be calculated. For con- 
venience of consideration, we will again turn to function 



(fc-i) 

Ak-l) _ 92. ..k 
J2...k — ■ 



(135) 



and using ( 35 1 we obtain the form suitable for the given 
case 



+00 



Jk-i) 

J2...k 



(136) 



ip{x) 



exp(— /3uoA;) — cxp(— /3uo) 
1 - exp{-l3vo) 



dx. 



Expanding the right side of ( 136 ) up to the terms of 
the second order with respect to /3u and making simple 
transformations we obtain 



+ 



i<j;i=2 



where 



and 



+00 



</>0i = / v{x)[va ~ v{x + Xi)]dx 



(137) 



(138) 



+ 00 



[v{x)vq — v{x + Xi)v{x + Xj)]dx. (139) 



Passing to u with the help of (34), we have 
z.(z,r)=/3^|^[(fc-l) / Z^^'i...fedr2...drfc (140) 

k=2 ' 

{k-l){k-2) f (fe) 

ha l^0 2...kdr2...drk 



Then, differentiating ( 140 1 with respect to the chemical 



potential and summing the series, we obtain a connection 
between the nonlinear part of the surface number density 
and the first localized distribution functions 



h2^o}dr2 + ]^ / (l>23 J^o^2,3dr2dr3 



(141) 



where we used representation (171, as well as relation ( 76 1 



It should be stressed that Ursell functions comprised in 



( 141 1 relate to the uniform medium. 



Let vo > 0, then it follows from ^ - (|27|) that for T ^ 
z' ^0 

6, ^ 9^ (142) 



In this situation, as one can see from (61 1, the coordi- 
nate of zero adsorption surface becomes 



u oo 

x, = - J eix)dx + J^ix)da 



while Henry constant of adsorption is 
KHix') = x' -xo{T), 



(143) 



(144) 



which follows from (|60|). For /i'^ J'\ 7 and Qs we obtaine 



from (|136|, (|118|), (|72f respectively, 

fe 





+ 00 




Ak-l) 
J2...k ^ 


j 6{x) 


1 - 




— 00 








00 


-i^^P{z){x' 


- Xo) + 


I 






-00 






-|-oo 


Qs = q{z){x 


- Xq) + 


/ 






— 00 



i-l[eix + x,) 



dx. 



(145) 



Expressions (1431 - (147), as expected, exactly cor- 



respond to the results for impermeable wall, obtained 
in 



20 



17. Static membranes 

Within the framework of the problem of permeable wall, 
the problem of static membrane can easily be solved. The 
static membrane will be understood as the region of time- 
independent field, which is a surface (that has, however, 
a certain thickness), with the field equal to zero on both 
sides of this region. In other words, the problem of per- 
meable wall will be considered for 



/3vo = 
I3v{x) ^ 0. 



(148) 



The basic expressions can be most easily obtained from 



the forms like (116). Assuming z' = z in (52 1, (116), we 
have 



and 



Qs^ J [g{x) - g] dx 

— 00 

1^ j [P-Pt{x)\dx, 



(149) 



(150) 



where g and P are the number density and the pressure 
of the fluid far from the membrane. 
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We obtain from ( 42 ) for this case 
O = -PV + "fA. 



(151) 



For the case of the uniform medium we have the system 
of equations ([9]), (10), and it may be written in the form 



In the first glance, it may be concluded from ( 151 1 that 



there is no uncertainly in separating the volume and sur- 
face terms for membranes. In reality, for "thick" mem- 
branes, as one can easily understand from the form of 7 



n=2 
00 

g = z + nbnz" 



n=2 



{ 150 1, this uncertainty is conserved: it is contained within -where 



7. Indeed, the "thick" membrane does not differ from the 
voluminous field considered above. 

Henry constant of adsorption takes on the form 



(156) 



(157) 



K 



H 



exp(— /3w) — lldx. 



(152) 



It should be stressed that integration in (157) is per- 
formed in infinite limits |24| . 
Excluding z we obtain 



k+l 



(158) 



which is most clearly seen from (57) 



k=l 



In agreement with logics, Kh is negative for positive 
potentials and vice versa. 



where 

iiemiaiH anu viutj versa. n 11 M 1 ■ 

. . . ^ „ _1 ^ 1 + Z^ i "^j)' TT (j 

If we start from the basic expression (71 ), we will have /S^, = ^jj — _[ J_ — 



for the first term in the right side of this expression 



00 

-[P{z)-P{z)\{x -xq) ^ QksT j [l-exp(-/3u)]dx. 

— 00 

(153) 

The term with x' becomes zero, which is logical: the 
disappearance of phases must mean the disappearance of 
interface. The value oi x^ co: in this case adsorption 
is always present. 

The nonlinear surface coefficient becomes in this case 

00 00 
v{z,z)= j [P~Pt{x)\dx~QkBT j [l-exp(-/3u)]da;, 

— 00 —CO 

(154) 

which can be obtained, for example, through expansion 



(159) 

are so-called irreducible cluster integrals. Summing is 
performed over all sets of nonnegative vrij meeting the 
requirement 

k+l 

Y,{3~l)m,^k (160) 

|4j p. 144], which provides the same dimensionality of 
the terms of sum (159). 



In the presence of the field, we must start from equa- 



tions ( 86 ) , ( 20 ) . Averaging them over the transition layer 



of (117) in series of Pvq. 



we obtain the system 

00 

pp^ = ze + ^dn{zeY 

n=2 

00 

Q = ze + ^ndn{z6Y' 



(161) 



Because equations (106), (107) remain true for this 



case, the expressions similar to the above-mentioned ones 
with the substitutions Pt — > P* and Pt Pst are valid 
too. 

At last, expression for the function /C^^^) takes the 
form 



where 



-(Tt 1 1 [\{^]^'C!ndxidr^...dr,,. (162) 



+ 00 



l2...k — 



1 - kip{x) - 9{x) Yl d{x + X,) dx, (155) 



Here integration over the coordinate of 1-th particle is 
made in the limits of the transition layer, while integration 
over the coordinates of other particles is performed in 
infinite limits. 



which follows from the expansion of ( 136 ) in series of (3vq 



Replacement of (155) into (34) gives (154) for the version 



Similarly to (33) we may integrate over xi. Making the 



with P* , if the order of integration is changed. 



18. Near-surface virial expansion 



change of variables x[ — xi, r ■ = — ri (i = 2, . . . n), 
using the invariance of with respect to translations 

and dropping the primes, we get 



!,("-l)7/(n) ^„ J„ 
"2...n "0,2...n"'^2---ar„, 



Equation ((Till) and the properties of (fT07|, allow us where 



to construct the near-surface analogs of virial expansion. 
The results of this building coincide with the results of [20] 
for an impermeable wall. 



n 

9{x)[Y[e{x + x^) 

i=2 



dx. 



(163) 



(164) 
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In view of the identity of the functional connection of 



P and Q in ( 156 ) on the one hand, and P* and g in ( 161 ) 



on the other hand, we can immediately write the analog 

(165) 



of (1581 



4. The general expression was obtained for specific sur- 
face f2-potential of the nonuniform system (7), it is 
composed of the parts that are linear and nonlinear 
with respect to pressure ( 71 1 . The linear part is de- 



PP* 



Q 



k=l 



termined by the product of pressure difference and 
the distance from the separating surface to the sur- 



where 



face of zero adsorption ( 59 ) , while the nonlinear part 

v) has several uni- 



(the nonlinear surface coefficient 



versal representations (33 1, (34), (90), (117) 



n 



(166) 

while requirement ( 160 ) is conserved in the previous form. 



5. The surface number density, similarly to 7, is sep- 
arated into the part linear with respect to density, 
which depends on the position of dividing surface. 



Finally, taking into account (107), (115), we may write 



and a universal nonlinear one ( 72 ) 



and 



fc=i 



k + l 



■k+l 



E 

fc=i 



-Ski 



■k+l 



(167) 



(168) 



6. Henry constant of adsorption depends on the position 
of the dividing surface and is determined by expres- 



sion (60) 



An expression for the nonlinear surface coefficient v 
was obtained in the form of the series in powers of the 
activity, which is an analog of the cluster expansion 
for pressure (34). The coefficients of the series are 



Expressions (167) and (168) are preferable than (165) 



both Pt and Pgt can be derived from computer experi- 
ments. In addition, P^ can be calculated because it may 



represented as the integrals of the products of Ursell 
factors and multipliers depending on the external po- 



tential (35) 



be represented in quadratures ( 108 ) 



From the viewpoint of the nontriviality of equalities. 



averaging intervals in (167), (168) are to be chosen min 



For macroscopically smooth field of arbitrary config- 
uration (section [8| the expression obtained for {79) 
agrees with the expression for the fi-potential of the 
nonuniform system (|80|). 



imal but such that the equalities (106), (107) would still 
be valid. 

Thus, a connection of Sk with is determined by equa- 9^ ^he approach through pressure tensor gives an ex- 



tion (166). For example, several first relations are 



pression for v in quadratures (117), which allows one 



51 - 2d2 

52 = 3(^3 — 6^2 
^3 = 4(^4 



to operate on 7 in whole ( 118 ) 



(169) 



„ . , , 80 
2M2d3 + -^4 



10. The approach through the cluster expansion gives the 
expression for v in quadratures when using the "local 
pressure" 



(90) 



11. 



and they are identical to the corresponding relations be- 
tween Pk and bj |4, p. 144]. 

19. Summary 



As an average over the transition region in the vicin- 
ity of the surface, the tangential component of the 



pressure tensor ( 108 ) coincides with the pressure act- 



ing on the transverse wall (91) and with the "local 
pressure" (86 ). This fact, as well as the arbitrary po- 



sition of the dividing surface, allows us to consider the 



1. 



The general expression was obtained for the fl- 
potential (grand potential) of the system into which 
a force field of a finite value was locally introduced 



new variants of Kirkwood-Buff formula (89), (102), 



(116) 



(32). This approach allows one to introduce an arbi- 



trary position of the separating surface ( 36 ) 



12. The il-potential of the system in the field of arbitary 
configuration is determined by the integral of "lo- 
cal pressure" over the volume ( 83 ) , ( 86 ) . In view of 



2. The volume and surface terms of the r2-potential can- 
not be separated unambiguously; their interpretation 
is determined by the shape of the interaction poten- 
tial between the body and the particles of the fluid 



13. 



(37) - (39) 



3. 



In particular, depending on the situation, the surface 
terms can start from a term which is either linear 



the previous item, it is also determined by analogous 
integrals of Pt and Pst, when these values may be 
introduced. 

As an average over the transition region in the vicin- 
ity of the surface, the relation between the number 
density and the tangential component of pressure 
tensor corresponds to the usual macroscopic relation 



or quadratic with respect to the activity (37), (39), 



between a number density and a pressure (114) 



which corresponds to the presence or absence of ad- 
sorption in the classical understanding of this term. 



14. Approaches through the cluster expansion and 
through the pressure tensor for the problems under 
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consideration are completely equivalent within the A. 2 Factors of partial localization 



domain of their existence (section 13 ) 



15. The solution of the problem of interface in the force 
field of a finite value possesses the symmetry with 
respect to the permutation of the field and the free 



As far as we know, these functions were introduced for 
the first time in 



26 



fluid (section 14) 



16. Within the problem under consideration, expressions 
for 7, KH^ Qs were obtained for static membranes 



Some particles involved in these functions do not cause 
the decay when moving away (the delocalized group) , and 
some of them - do (the localized one). 

Introduce the notation 



(section 17 1 



„{m,k) 
'-'l...m+fc' 



(A.3) 



17. Within the problem under consideration, the para- 
metric relation between Henry constants of adsorp- 
tion and absorption was obtained (section [6]). 

18. At high temperatures, the nonlinear part of the sur- 
face number density is expressed through the func- 
tionals of Ursell functions of the 2-nd and 3-rd order 



where the superscripts m and k define the quantity of 
delocalized and localized particles, respectively (m — 
1, 2, 3, . . . , fc = 0, 1, 2, . . . ). The subscripts denote the co- 
ordinates of particles, with the first m particles being con- 
sidered delocalized, and the rest - localized. 

These functions are similar in structure to the Ursell 
factors of fc -t- 1-th rank, with the proviso that in the con- 



(1411 



struction by the type of (A.ll the first m particles (de- 
localized) are treated as a single compound particle. In 

(m,A:) 



19. The near-surface virial expansion (section [Ts]) gives other words, define '^-^j. by the equality 



the equation of state for the "two-dimensional" fluid 
in the region near the boundary. 

Appendix A Factors 
A.l Ursell factors 

Ursell factors can be defined by equality: 

I 

{n} a=l 



1 ^ < fc, ^ ^ fco 



fc, exp{-l3U^) = 1, (A.l) 



where {n} designates some partition of the given set of fc 
particles with coordinates ri, . . . , rj, into disjoint groups 
{no,}, / is the number of groups of given partition, ka is 
the size of the group with the number a. Summing is 
made over all the possible partitions; the sense of the re- 
quirement exp(— /3C/^) = 1 is evident: groups containing a 
sole particle do not give any contribution into the product 
in this case. 

(k) 



For example, some initial \ are 



U['^ =1 

= eM-PUl2.z) - exp(-/3C/L) 
- exp(-/3;7i%) - exp(-/?C/|,3) + 2 



(A.2) 



It is well known that one of the generating functions 



for Ursell factors is logarithm 22 
for the case of fc 



However, dividing the 
1 by Sy in the form of (|6]), 



series m 

one can easily make sure that there is another generating 
function - a fractional function, in the form of the ratio 
of series (Isl). This follows from the well known expansion 



for the number density of the uniform system ( 10 ), which 
is true far from the boundaries of the system. 

The existence of the fractional generating function will 
be proved in Appendix |B.1[ 



{n} 



X nexp(-/3C/'^'°+(™-^)'°-({na})), (A.4) 

a=l 

I 

l<fca<fc + l; ^ka = k + l; exp(-^f/^) = 1, 

a=l 



where designations are analogous to (A.l) on condition 



that summation is taken over all possible partitions of 
the set of fc -f 1 particles among which one particle is com- 
pound. 5aT is the Kronecker delta, and r is the number 
of the group involving a compound particle. 

is the distribution 



2{m.k) 



The generating function for B\ 

function Qq^\ ^ of GCE type. Expanding the partition 
function in ([8| and dividing the series, we obtain ( 13 ). 
The proof of this relation is given in Appendix |B.1[ 



The first in localized group are 



'^l...rn 

"l...m,-|-l 

„(m,2) 

'~'l...m+2 



= exp(-^C/™.„J 
= exp(-/3C/™+i+ 
= exp(-/3C/™+^^^ 

X e^v{-PU'^+l,ra+2) 



exp(-/3C/{"..„) 
exp(-/3C/{".+^+i) 

+ 2exp(-/3C/™.™) 



(A.5) 



and in delocalized one are 

including, for homogeneous medium 

= 1. 



(A.6) 



(A.7) 



It is obvious from (A.5), and (A.6) that factors of par- 



tial localization generalize the concepts of Boltzmann fac- 
tors and Ursell factors, including them as the limiting 
cases. 



For factors B\ 



(m^k) 



.m+k 



a number of recurrence relations 



hold, which ensure the existence of various physical links; 
some of them are given in Appendix [B| 
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Appendix B 



Recurrence relations for 

{m,k) 
m+k 



Many relations and operations of statistical mechanics 
are provided by definite classes of recurrence relations for 
^i^'m+k- ^'^^ brevity, we shall say that an operation gen- 
erates a recurrence relation or a class. In this contribution 
only some of them will be considered. More detailed in- 



formation can be found in 27 



Appendix C Proof of (Pt) = (P*) 

Let us prove the equality of the integrals of P* and Pt 
taken over the entire volume. To prove this equality, con- 
sider the system inside the volume shaped as rectangu- 
lar parallelepiped strongly stretched along the y axis (see 
Figj2]). The interface is inside the system. The signifi- 
cance of field/fluid interaction is still defined by the func- 
tions 9. 



B.l Correspondence to the definition 



Equating expressions ([8|) and (13 1, expanding the series 
for Sy, and performing multiplication of the series, we 
arrive at 



"l...m+fc ~ "l...m+fc 



EE 

n—l samp 



n+1. ..m-f A; ' 



(B.l) 




where m > 1, fc > 1, and the internal sum is taken over 



samplings n of initial localized k particles. Relation (B.l ) 



Figure 2: 



IP* 



is proved either by direct exhaustion of partitions in ac- 



cordance with ([XIl), or by repeatedly substitutions the Let us integrate expressions (|86| and ( |108[ ) with respect 



expression for in the right-hand side of (B.l I. 

Paper 22 gives another recurrence relation for Ursell 



functions which in terms of B 



l...m+k 



looks like 



(i,k) 
i...fc+i 



B 



(A;+1,0) 
l...k+l 



k 

E 

n=l 



K(">0)ra(l,fe- 



.k+l 



perm 

(B.2) 

where fc > 1, and square brackets denote averaging over 
permutations of particles. 

1 



to the coordinates of the first particle. The integration 
is performed over the whole system volume denoted by 
dashed line in Fig. [2j We put that this volume is defined 
by the characteristic function ■01. 

Expand the integrands into a series of the activity using 
([22|). Equating the coefficients at equal powers of z, we 



have 



In principle, this is the same equation (B.l) at ra 



» fe+2 



(C.l) 



but written in a symmetric form. Note that in the given 
case there is no necessity in symmetrization - equation 
(B.l ) holds rigorously in asymmetric form as well. Never- 



(fc + 2) 



. fc+2 



^1 [ n (2/2 - yi)B^i.:''k+2^^dri...drk+2, 



i=l 



theless, (B.l) can also be brought into a symmetric form. 



however, this calls for averaging over permutations of both 
the sums and the left-hand side of the equation. 

B.2 BBGKI equation 

In the presence of external field this set of linking equa- 
tions may be written as j4| p. 205] 



where fc = 0, 1, ... . Here we equated the coeffi cients at 
The kind of the recurrence relation (B.4) required 



to prove (C.l) is of the form 



fc-t-2 



(C.2) 



i=2 

-P j el^.m+iViUim+idrm+i. (B.3) 



This relation, after substituting expressions (22) in it 



and equating coefficients at equal powers of z, generates 
the differential recurrence relation 



Multiplying \G.2\ by (?/2 - yi)tpidi ■ ■■0k+2 and inte- 
grating over all coordinates, we obtain using the symme- 
try under permutations inside localized and delocalized 
groups of ;B('"'") 

-keT 9^ (2/2 - yi)^M'.l+2 '^^i-'^^fe+2 

•' i=i ^2 

^h + kh + kh, (C.3) 



i=2 



rn-\-n 



E Mr.rC'^Vr^H, (B.4) 



i— m+1 



which we use to prove (107), (113). In the second term 



where 

h - 

h = 



_ fc+2 



i^i[X{9^{y2-yi)Bt{ 



fc) du2i , , 
k+2- — dri...drk+2, 



on the right all permutations of particles of the localized 
group (n) with the second particle are exhausted. 



i=i 

_ k+2 



dy2 



(C.4) 



^1 [ n (2^2 - y3)B2x..k+2^^^dri...drk+2, 



i=l 



(C.5) 
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k+2 



k) 



du2l 



— dri...drk+2- 

+ dy2 

(C.6) 

Here in the integrals I2 and the change of variables 
ri O rs is performed. 



The integral in the left-hand side of equality ( C.3 ) may 



be taken by parts, since here factors 9i depend solely on 
Xi. Wherein, due to locality of Ursell factors, and factors 
g(i,n) ^j.g exactly the factors of this kind, the integral 
term goes to zero. Thus, 

i'l [ n - yi)Qy-^^iAX'2dri-drk+2 

^ k+2 

J B[':'^l'^dr,...dr,+2. (C.7) 



It is readily seen that 

For this purpose, calculate the difference 

. k+2 



(C.8) 



V'l [ n ^] iy^-yi)^fx..k+2^^d'ri-drk+2 



i=l 



dy2 



(C.9) 

Making the change of variables ri ++ r2, we see that 



(C 



the right-hand side of (C.9 1 coincides with I2, so we have 



Finally, consider the integral I3. Obviously, the factor 
("03 — ipi) is nonzero only in the case when the first and 
the third particles are on different sides of the system 
boundary. It is seen better from the expression 



(C.IO) 



where we introduce the function Xi = ^ ~ ''Pi- 

When the 1-st or the 3-rd particle moves away from 
the boundary of the system, the integrand decays rapidly 



owing to local character of the product B^2i'\+2'^''^'^i/'^y'^- 
Hence, the integral is proportional to the surface of the 
parallelepiped at hand, but not all of its faces make a 
contribution to /s. 

Consider the top face of the system, where due to re- 
moteness from the interface we can consider that 0i — 
1,1 = 1,2,3.... Making the substitution (C.IO I, we ar- 



rive at the conclusion that condition = demands in 
this case that 



^3Xiiy2 - y3)B^2.i'lk 



+2 



du2i 
dy2 



dri...drk^ 

du2i 



(C.ll) 



= / ^iX3{y2-y3)B^2f^^k+2^^dri...drk+2- 



However, obviously, equation (C.ll) is true by symme- 



try when the y axis is parallel to the plane defined by 
the function ip. To prove this it suffices to consider the 
mirror reflection operation with its mirror-plane coincid- 
ing with the upper face of the system. Thus, we may 
conclude that contribution to I3 from the whole region 
of homogeneity above the interface equal zero except the 
side faces y — const. In the region of homogeneity below 
the interface situation is completely analogous due to the 
fact that in this region 9i = const, i — 1, 2, 3 ... . 

So, only the side faces y = const and two narrow strips 
on the faces z = const make contributions to I^. The 
contribution of these strips is of higher order of smallness 
(proportional the linear sizes of the system) and may be 
neglected. Thus, we have the estimation 



A 
s 



0{Szy). 



(C.12) 



The contribution of I3 tends to zero with the system 
length (along the y axis) tending to infinity. Here Szy is 
the area of the face zy. 



In view of (C.8I and (C.12), equation (C.3) gives (C.I). 
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